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Abstract: Employing the quantum Hamiltonian describing the interaction of two-mode light (signal- 
idler modes) generated by a nondegenerate parametric oscillator (NDPO) with two uncorrelated squeezed 
vacuum reservoirs (USVR), we derive the master equation. The corresponding Fokkcr-Planck equation 
for the Q-function is then solved employing a propagator method developed by K. Fesseha [J. Math. 
Phys. 33 2179(1992)]. Making use of this Q-function, we calculate the quadrature fluctuations of the 
optical system. From these results we infer that the signal-idler modes are in squeezed states. When the 
NDPO operates below threshold we show that, for a large squeezing parameter, a squeezing amounting 
to a noise suppression approaching 100% below the vacuum level in the of the quadratures can be achieved. 



1 Introduction 

Nonclassical effects of light such as squeezing, antibunching and sub-Poissonian statistics have been 
attracting the attention of several authors in quantum optics over the last decades [1, 2, 3, 4, 5, 6, 7, 8, 9]. 
A review article on nonclassical states of the first 75 years is found in Ref. [2] . 

Squeezed states are nonclassical states characterised by a reduction of quantum fluctuations (noise) 
in one quadrature component below the vacuum level (quantum standard limit), or below that achievable 
in a coherent state [4, 6] at the expense of increased fluctuations in the other component such that the 
product of these fluctuations still obeys the uncertainty relation [4, 5]. 

It was Takahashi [10] who, in 1965, first pointed out that a degenerate parametric amplifier enhances 
the noise in one quadrature component and attenuates it in the other quadrature. This prediction has 
been confirmed by several authors for degenerate and nondegenerate parametric amplifiers and oscillators. 
Operating below threshold, the parametric amplifier is a source of squeezed states. In the initial experi- 
ments carried out to observe squeezing, a noise reduction of 4-17% relative to the quantum standard limit 
has been obtained [11]. In order to increase the gain, the parametric medium may be placed inside an 
optical cavity where it is coherently pumped and becomes a parametric oscillator [12, 13, 14, 15, 16, 17]. 

An optical parametric oscillator is a quantum device with a definite threshold for self sustained 
oscillations. It is one of the most interesting and well characterised optical devices in quantum optics. 
This simple dissipative quantum system plays an important role in the study of squeezed states. In a 
parametric oscillator a strong pump photon interacts with a nonlinear-medium (crystal) inside a cavity 
and is down-converted into two photons of smaller frequencies. In the NDPO we assume that the strong 
pump photon is down converted into two modes and these modes are referred to as signal and idler modes. 

A quantum-mechanical treatment of different optical systems such as the NDPO is essential as they 
may generate squeezed states with nonclassical properties which have potential applications in quantum 
optical communications [4] and computation [19], gravitational wave detection [19, 20, 21, 22], interfer- 
ometry [22, 23, 24], spectroscopical measurements [25] and for the study of fundamental concepts. 

For systems with nonclassical features such as the NDPO, for which the Glauber-Sudarshan P-function 
is highly singular [8, 26, 27], one may use the Q-function. The Q-function is expressible in terms of the 
Q-function propagator and the initial Q-function. It is possible to determine the Q-function propagator 



by directly solving the Fokker-Planck equation. In this paper, we find it convenient to evaluate the 
Q-function propagator applying the method developed in [28] . 

The main aim of this paper is to calculate the amount of squeezing that could be generated by the 
NDPO coupled to two USVR with the help of the Q-function. We show that squeezing of the output 
may be optimised and reach 100%. 



2 The Master Equation 

The description of system-reservoir interactions via the master equation is a standard technique in quan- 
tum optics [6, 7]. In this section, however, we found it useful to include a non-detailed a derivation of 
the master equation describing the interaction of the signal-idler modes generated by a NDPO coupled 
to two USVR in order to make the paper more self-contained. 

Denoting the density operator of the optical system and the squeezed reservoir modes by x(t), the 
density operator for the system alone is defined by 

p(t) = Tr R (x(t)), 

where Tr^ indicates that the trace is taken over the reservoir variables only. The density operator x(t) 
evolves in time according to 



dt 



(1) 



where Hs R (t) is the Hamiltonian describing the interaction between the system and the reservoirs. Note 
that the Hamiltonian describing only the reservoir modes (H R (t)) is not involved in the derivation as 
it cancels out when we apply the cyclic property of the trace. Furthermore, in oder to simplify our 
calculations, the Hamiltonian that describes the interaction of the system with the pump mode (Hs{t)) 
will be added at the end of the derivation. 

Since initially the system and the reservoirs are uncorrelated, one can write, for the density operator 
of the system and the reservoirs at the initial time (t = 0), that x(0) = p(0) <8> R [6], where /3(0) and R 
are the density operators of the system and the reservoirs at the initial time, respectively. Then in view 
of this relation, Eq. (1) results in 



dm 

dt 



1 

ih 



H SR {t),p{0)®R 



1 



dt' 



H S R(t),[H SR (t'), x (t')] 



(2) 



Applying the weak coupling approximation which implies that x(i') = p(t') ® R [6], it follows that 



df>(t) = ]_ Tir 
dt ih 



H SR (t),p{0)®R 



- I dt'Tr, 



H SR (t),[H SR (t r ),p(t')®R] 



(3) 



We consider the system to be a two-mode light with frequencies uj a and u>b in a cavity coupled to two 
USVR. The interaction between the two-mode light and the squeezed vacuum reservoirs can be described, 
in the interaction picture, by the Hamiltonian 

H SR (t) =ih\j[^\j{tfA j e i ^- w ^ t -aA]e- i ^-^^^ , (4) 



in which a (d^) and b (&) are the annihilation (creation) operators for the intracavity modes and Aj {Aj) 

and Bk {R\) are ^ ne annihilation (creation) operators for the reservoir modes with frequencies ojj and 
u>ki respectively. The coefficients Aj and Afe are coupling constants describing the interaction between 
the intracavity modes and the reservoir modes. Applying the cyclic property of trace and the relation 
Tr R (R (g> Hsn(t)) = (HsR(t)) R , and taking into account that, for squeezed vacuum reservoirs [6], 

(A j ) R = (Al) R = (B k ) R = (Bl) R = 0, 



one can show 



ih 



-TV, 



H SR (t),p(0)®R 



0. 
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and as a result, expression (4) reduces to 



dm 

dt 



1 



^ / dt'Tv R { H SR (t),[H S R(t'),p(t')®R] 



(5) 



Applying the Markov approximation, in which p(t') is replaced by p(t), and using the cyclic property of 
the trace, the above equation can be expressed as 



dm 

dt 



1 



dt' 



(H SR (t)H SR (t')) R p(t) - (H SR (t')H SR (t)) R p(t) 



~p(t)(H S R(t)H SR (t r )) R + p(t)(H S R(t')H SR (t)) R 
We note again that for squeezed vacuum reservoirs [6] 

(AjA^h = -M A Si, 2 j a -j, 
(itii) fl - N A S jtl , 
(A j Al) R = (N A + l)S jtl , 

where Sj.i is the Kroneckcr delta symbol and 



{AjB^r = (AjBlU = (AW m ) R = 0. 



(0) 

(7a) 
(7b) 
(7c) 

(8) 



This equation is a consequence of the fact that the two squeezed vacuum reservoirs are uncorrelated. 
The parameters N A , N R , M A and M R describe the effects of squeezing of the reservoir modes. Actually, 
the parameters N and M represent the mean photon number and the phase property of the reservoirs, 
respectively, and are related as |M| 2 = N(N + 1). Furthermore, introducing the density of states g(u), 
where 



Jo 



dw g{u})\{uj)X(2oj a — w), 



and setting t — t' = r, one can easily show that 



/' 

Jo 



dt' e ±i ( w "- u )(*- t ') 



dr e ±%[u 



(9) 



Since the exponential is a rapidly decaying function of time, the upper limit of integration can be extended 
to infinity. Making use of the approximate relation 



=±iC«-.-")T =7r< j( Wo _ w ) j 



and applying the property of the Dirac delta function to the integrals of Eq. (6), we get 



irg(uj a )\ 2 (uj a ) 



2 ' 



(10) 



(11) 



where j A = 2Trg(uj a )X 2 (ui a ) is the cavity damping constant for mode A. Similarly one can also show that 
the cavity damping constant for mode B is given by 75 = 2itg(uii l )\ 2 (uib). 

In view of Eqs. (7-11), after evaluating lengthy but straightforward consecutive integrations, Eq. (6) 
takes the form 



dm 

dt 



= -^-(N A + 1) 2a / 5(i)a t - tfap{t) ~ p{t)a)a 



1 A N A 



2d) p(t)a — da) p(t) — p{t)dd) 



1A.M A 



2a t / 3(t)a t + 2ap{t)a - a t2 /5(i) - p(t)a n - a 2 p(t) - /5(t)a 2 



IB 

' 2 

jbM b 



+ 1 f(N B + 1) [2 bp{t)tt - b r bp{t) - p{0b 



1bN b 



2# mi - itf p{t) ~ m)t> h] 



2 pp(tp + 2 im h - b t2 m - mtf 2 - b 2 m - P{t)b 2 



(12) 
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In the cavity we consider two-modes of light known as the signal and idler modes produced by the 
NDPO. The cavity has one single-port mirror in which light can enter or leave through while its other 
side is a mirror through which light may enter but can not leave. In this system we assume that a strong 
pump light of frequency c^o interacts with a nonlinear-medium (crystal) inside the cavity and gives rise 
to a two- mode squeezed light (the signal-idler modes) with frequencies tv a and oj^ such that ujq = uj a + uj(, . 
With the pump mode treated classically (the amplitude of the pump mode is assumed to be real and 
constant), the interaction of the system is described, in the interaction picture, by the Hamiltonian 



(13) 



where k is the coupling constant and 70 is the amplitude of the pump mode. Hence the master equation 
for the NDPO coupled to USVR, in view of Eqs. (12) and (13), takes the form 



dt 



abp(t) - p(t)ab + p(t)aW - a ] P p(t) 



JA 
' 2 
lA M A 



+-^-(N A + 1) 2a / 5(i)a t - a)ap(t) - p(t)tfa 



1aN a 



2 a) p(t)a — aa) p(t) — p(t)aa) 



2 tfp{t)tf + 2 ap(t)a - a) 2 p{t) - p{t)a) 2 - a 2 p(t) - p(t)a 2 



IB 

' 2 
1bM b 



(N b + 1) 2bp(tp - &bp(t) - p(t)b r b 



1bN b 



2tfp(t)b-btfp(t)-p(t)bP 



2 tfp{tp + 2 bp(t)b - tf 2 p(t) - p(tp 2 - b 2 p(t) - p(t)b 2 



(14) 



This equation is the basis of our analysis and describes the interactions inside the cavity as well the 
interaction of the signal-idler modes produced by the NDPO and the squeezed vacuum reservoirs via the 
partially transmitting mirror. This master equation is consistent to that given in Ref. [7] except that the 
expression there is for a single mode in a cavity coupled to a single mode vacuum reservoir. 



3 The Fokker-Planck Equation 



In this section we derive the Fokker-Planck equation for the Q-function. In order to obtain the Fokker- 
Planck equation for the Q-function corresponding to the master equation (14), one has first to put all 
terms in normal order. Applying the commutation relations 



a, /(a, a*) 
a\f(d,a)) 



Sat ' 
d/(a,qt) 
da 



(15a) 
(15b) 



one can verify that ap = pa- 



cv p- 



dp 



|, where the density operator p = p(a, a,*, t) is considered 



to be in normal order. Making use of Eqs. (15), the relation [a, be] = b[a,c] + [a,b]c and the photonic 



commutation relation ad 
dp(t) 



+ 1, the master equation (14) can be written as 



dt 



-«7o 

1A 



dp . . ^dp 
—z-a + a' — 
96t db 

d 



St 



2 (tfx + UI^) 



dp dp t d 2 p 

— H —b J — 

da da* 



d 
<9at 



d 2 p 



(atp) 



dadb 
n d 2 P 



~~2~ Ma 



&p 
da 2 



d 2 p 



<9dt2 



IB 

' 2 

IB 



{N B + 1) 
d 2 p 



\d_ 

Ob 



dtf 



ftp) 



dada^ 



d 2 p 



datdtf 
1aN a 



^p) 



dii 



(pa) 



dbdW 



JbN b 



ft P ) + - i (pb) 
db 



iB 



d 2 p 



db 2 9at2 



(16) 
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In order to transform this equation into a c-numbcr Fokkcr-Planck equation for the Q-function, one needs 
to multiply it on the left by (a, (3 | and on the right by | a, P), so that 



dQ 
dt 



i 92 
«7o ( ^ n „ + 



if 







dadp da* dp* dp* 
d 2 , 1A / d 
2 y da 
1b , d 



+1a{N a + 1) 
+1b{N b + 1) 



dada* 
d 2 



d . 
dp a 
d 
da* 
d 



d ^ d , 
da* da 
s 1A M A , d 2 



d 2 



'■da 2 da* 2 ^ 



dpdp* 



where 



2 K dp l3+ dp* 13 ' 



1 



7s M b , 9 



d 2 



dp 2 dp* 



Q, 



(17) 



Q = Q(a*,a,/3*,/3,t) = —~{a,P \ p(a\a,b\b,t) | a,/?). 



Expression (17) is the Fokker-Planck equation for the Q-function for the signal-idler modes produced by 
the NDPO coupled to two USVR. To obtain the solution of this equation, we introduce the Cartesian 
coordinates defined by 



a = x\ + iyi, a 



and note that 



One can show that 



X\ - 


- iyi 


P = 


X2 + iy2, P* = x 2 




(a- 


-a ) 


X2 = 






d 




, d 


■ 9 \ 




da 




"dx\ 


dy\ 




d 




, d 


■ 9 \ 




dp 


-\ 


"dx 2 


dy2 





(18) 
(19) 

(20a) 
(20b) 



Thus combining these results and their complex conjugates, one readily obtains 



dQ 
dt 



2 



d 2 



2 \dxidx 2 dyidy 2 

1a {n a + 1) / d 2 d 2 



«7o 



d d 

-X2 + 7: X\ 



4 

1b{N b 



1) 



dx\ dy\ 
d 2 d 2 



dy 2 2 



dx\ 



dx\ dx2 
jaMa ( d 2 d 2 
4 \dx 2 dy\ 

JbMb ( d 2 d 2 
4 



\dx 2 2 



dy\ 



d d 
dyi dy 2 



Q, 



(21) 



where Q = Q{xi, x 2 , yi, 2/2, *)• 

Next, introducing the transformation defined by x\ 
can verify that 

x = 2^ Xl +' T2 )' S 

y = -Avi ^ 2/2), 

In view of these relations, it follows that 



x + u, X2 = x — u, yi = y + v, y 2 = v — y, one 



u = -{x\ - x 2 ), 



d 


1 


- d 


d 


d 


1 


- d 


d ■ 


dx\ 


~~ 2 


-dx 


du 


dx 2 


~~ 2 


-dx 


du. 


d 


1 


- d 


d ■ 


d 


1 


d 


d ■ 


dyi 


~~ 2 


.dy 


dv. 


dy2 


~~ 2 


dv 


dy. 



Making use of Eqs. (22, 23) in Eq. (21) and setting ~/ A = j B = 7, N A = N B = N and M A 
for convenience, one arrives at 



Mf 



(22a) 
(22b) 

(23a) 

(23b) 
= M, 



dQ 
dt 



«7o + 7(iV - M + 1) d 2 k 7o + j{N + M + 1) d 2 k~/ - j{N - M + 1) d 2 



8 dx 2 
K7o ~ l(N + M + 1) d 2 2k 7o + 7,9 



dy 2 



d 



8 

2k7o -7/9 



d 



8 dv 2 2 dx dy 2 y du ' dv 

which is the Fokker-Planck equation for the Q-function where Q = Q(x, y, u, v, t). 



du 2 

Q, (24) 
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4 Solution of the Fokker-Planck Equation 



In this section the explicit expression for the Q-function that describes the optical system is derived. In 
order to solve the differential equation (24) using the propagator method discussed in Ref. [1] , one needs 
to transform the above equation into a Schrodinger-type equation. This can be achieved upon replacing 
(lk'Tk'S^>Tk> x >y> u > v ) &ndQ(x,y,u,v,t) by (ip x ,ip y ,ip u ,ip v ,x,y,u,v) and | Q(t)) respectively. Hence 
Eq. (24) can be expressed as 



: d\Q(t)) 
' dt 



A 



1 ~2 A 2 ~2 



A3 -2 



H I Q(t)), 



A4 ~ 2 , • ^5 , A „ 



■PyV) -*y(Pu" 



P*v) \Q(t)) 



(25) 



where 



Ai, 2 = k 7o + 7 (7VtA^+1): 
A 3 .4 = K 7o - 7 (A r TA^+l): 
A 5 . 6 = 2k 7o ± 7. 



(26a) 
(26b) 
(26c) 



A formal solution of Eq. (25) can be put in the form 

I Q(t))=u(t) I Q(0)>, 
where u(i) = exp(—iHt/hbar) is a unitary operator and 

. A2 .0 . As _o . A4 A5 



H 



~ + + _ -7T\P* X +Pyy) + TriPuU + Pvv) 



Ag 

8 rv 2 v "~ ' r> "" ' 2 
is a quadratic quantum Hamiltonian. Multiplying (27) by (x,y,u,V | on the left yields 

Q(x,y,u,v,t) = (x,y,u,v | u(t) \ Q(0)), 



(27) 



(28) 



(29) 



where Q(x,y,u,v,t) = (x, y, u, v \ Q(t)}. Introducing a four-dimensional completeness relation for the 
position eigenstates I = J dx' dy' du' dv' \ x' , y' , u' , v') (x' , y' , u' , v' | in expression (29), one can see that 



where 



Q(x, y, u, v , t) = I dx dy du dv Q(x, y, u, v, t\x ,y ,u ,v , 0)Q o (x ,y ,u ,v) 



Q (x',y',u',v') = {x',y',u',v'\Q(o)) 



(30) 



(31) 



is the initial Q-function and Q(x, y, u, v, t\x', y', u', v', 0) = (x, y, u, v\u(t)\x', y' , u', v') is the Q-function 
propagator. 

Following Fesseha [1], the propagator associated with a quadratic Hamiltonian of the form 



H(xx, ...,X n ,pi, ...,p n ,t) = ^ [ a tPi + h{t)f>iXi + c t (t)x 2 

i=i 



is expressible as 

Q(x\ 1 . . ., x n , t\xi , .. ., x n , 0) 



jj 2 
2^ 



n 



d 2 S r _ 



dxjdx'j 



T ex P 



/ bj{t')dt' + iS c 



(32) 



(33) 



where S c is the classical action, £ is a parameter related with operator ordering and is constant 
different from zero for the Hamiltonian to remain quadratic. Comparing Eqs. (32) and (28), it follows 



that a,i = 



jAi,2,3,4, (xi, X2, X3,Xi) = (x,y, U, v), C x = C y 



-*u — L i' 



0, b x 



A4 

2 ' 



b u = b v = 4^ and the antistandard operator ordering £, = \- Thus the Q-function propagator associated 
with the Hamiltonian (28) is expressible as 



Q(x,y,u,v,t\x',y',u',v',Q) 



4tt 2 



d 2 S r d 2 S r d 2 S r d 2 S r 



dx'dx dy'dy duidu dv'dv 



»s c + (A5 ; Ae) t 



(34) 



(i 



In order to obtain the explicit form of this expression, one has first to determine the classical action. To 
this end. the Hamiltonian function corresponding to the quantum Hamiltonian (28) is given by 



H 



2 2 , -^3 2 



. A4 2 A5 



A fi 



' l Y Px ~ l ~8 Py ~8 Pu ~8 Pv ~ PyV ' ~2^ Pu P 



(35) 



With the help of the Lagrangian L = ^2i i XiPi — H and the Hamilton equations ii = [i = x,y,u,v) 
one can readily show that 



2i , . A5 ,2 , 2i , A5 ,2 2i , . A 6 ^2 2i,. A 6 ^2 

l = y 1 \ x + — x ) +—yy + —y) -—( u -— u ) -—( v -— v )- 



2 ' A 2 
Applying the Eulcr-Lagrange equations 



along with Eq. (36), leads to 



d f dL\ dL 
dt V dx,; ) dxi 



-A 5 > 2 



= 0, 



(36) 



(37) 



A5 \2 



-A 6n 2 



A 6x 2 



The solutions of these differential equations can be written as 

x(t) = a\e 2 + a 2 e 2 , 
u(f) = cie 2 + c 2 e 2 % 



v(t) = die—* + d 2 e~~. 



(38a) 
(38b) 



Now substituting these expressions and their corresponding first order time derivatives into Eq. (36), 
the Lagrangian takes the form 



Ai A2 



A3 Xa 



-A 6 i 



On account of the above result, the classical action defined by S c = fl L(t)dt takes the form 



(39) 



Applying the boundary conditions Xi(0) = x\ and £j(T) = x" in Eq. (38), one can obtain that 



ai 



II CUT 1 

x e 2 — x 



61 



11 ±S-T 1 

V e 2 - y' 



,A R T 



C2 = 



we 2 



d 2 = 



e -KT_ 1 > " e- A ^ T -f ' 
Inserting the above expressions into Eq. (39) and replacing (x", y", u" , v" , T) by (x,y,u,V,t) yields 



S r = 2iXc 



(x' — e i*) 2 (y' — e 2 '*) 



A 5 t\ 2 



Ai(e A =*-l) A 2 (e A =*-l) 



2iX f 



[u — e 2 



' A 6 f x2 



w — e 2 



- A 6 f s 2 



A 3 ( e -A6t_i) A 4 (e- A «*-l) 



and employing this relation the following results are obtained: 



d 2 s c 

dxdx' 

d 2 S c 
dudu' 



4zA 5 e^* 
Ai(e^*-1)' 



4iX 



8 2 S C 
dydy' 

d 2 S c 
dvdv 1 



MXseJ^ 
~X 2 {eM-iy 

MX§e~~$~ t 
A 4 (e~ A ^ - f) 



(40) 

(41a) 
(41b) 
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Thus, in view of Eq. (41). the Q-function propagator (34) takes the form 

4A 5 A 6 e (A5-Ae)t 



Q(x,y,u,v,t\x ,y ,u ,v ,0) 



7rVAiA 2 A 3 A 4 (eM - l)( e -*«* - 1) 

2A 5 fx 12 -2xe b T t x l + x 2 e x '' t 



x cxp 



(e^* - 1) V 
y' 2 - 2ye^ Lt y' + y 2 e X5t 



Ai 



2A fi 



( e -A 6 t _ X ) 



u' 2 - 2ue~^ l v! + u 2 e~~ Xet 



*„.' i „,2„-A 6 t 



v 2 — 2v e 2 v -|- 
A4 



(42) 



Considering the signal-idler modes produced by the NDPO to be initially in a two-mode vacuum 
state, the initial Q-function is expressible as 

Qo(a',/9') = \{a',/3'\0,0)(0,0\a',p') = exp(-a'V - 
and in terms of the Cartesian variables of expression (19), this equation becomes 
Qo(4» 4, V'l, y'2) = \ exp [ - (a£ + x'i + y' 2 + y 12 ) 



Furthermore, in terms of x' , ?/, m' and v' , one can write 



dx x dx 2 dy 2 Qo(£i, z 2 , 2/u 2/ 2 ) = / dx d V du dv Qo(x ,y ,u,v), 



where 



Q (x',y',u',v') = -^r exp 



2 (a; + y + w +wj 



and J is the Jacobian of the transformation of x\, x 2 , y\ and j/ 2 with respect to x, y, y and v. Making 
use of Eq. (19) in the Jacobian, one can show that \ J\ = 4. Hence 



Qo(x', y', u', «') = — exp - 2(x' 2 + y' 2 + u' 2 + v' 2 ) 



(43) 



Substituting expression (41) into Eq. (34) and then combining the result with Eq. (43) and finally carrying 
out the integration in Eq. (30) applying the relation 



dx 1 exp 



kx' 2 + dx' 



4A: 



k > 0, 



the Q-function takes the compact form 

4 



Q(x,y,u,v,t) 



7T v^i^aiol 



cxp 



— a; 2 y 2 

di a 2 a 3 a 4 



where 



ai,2 = 



A12 e A5t -l +A 5 



Ct3,4 = 



A 3 , 4 (e- A6 * - 1) + A 6 



A 6 e _A6t 



-u 2 - -v 2 



(44) 

(45a) 
(45b) 



It can be easily verified that the Jacobian of the inverse transformation is | J'| = |. One can then write 
J dxdy dudvQ{x, y,u, v, t) = I dxi dx 2 dyi dyi Q'(xi,X2,yi,y2,t), 
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in which the final expression for Q'(x\, x%, Hi, JJ2, t) is obtained from Eq. (44) employing the inverse 
transformations (22). Upon carrying out further inverse transformations (19), the required final form 
the Q-function for the signal-idler modes produced by the nondegenerate parametric oscillator (NDPO) 
coupled to two uncorrclated squeezed vacuum reservoirs takes the form 



Q(a,a*,/3,l3*,t) = ^ exp [ - 6i(|«| 2 + |/3| 2 ) + b 2 (a(3 + a* f3*) + b 3 (a/3* + a*/3) 
+ b -±(a 2 + a* 2 + P 2 + (3* 2 ) 



where 



D = 

61,2 = 
63,4 = 



1 



^aia 2 a 3 a 4 
1 " 



1 

±ai 



± 



1 


1 

\ 


1 " 


0,2 


03 


04 


1 


±1 


1 






T — 


0-2 


«3 


04 



(46) 

(47a) 
(47b) 
(47c) 



This Q-function is useful to calculate the expectation values of antinormally ordered operators and 
consequently the quadrature variances. It could also be used to calculate the photon number distribution 
of different optical systems. In this paper, this function is used to calculate the quadrature fluctuations 
(variances) of the NDPO coupled to two USVR. It can be readily verified that the Q-function (46) is 
positive and normalised. 

Now we proceed to obtain the expressions for the Q-function for some special cases of interest: For 
the case when there are no squeezed vacuum reservoirs (r = 0), that is, when the external environment 
is an ordinary vacuum, the Q-function (46) takes the form 



Q(a, a*, /?,/?*,*) 



1 



7r z ai<23 



■ exp 



^)<'»' 2 + '<>' 2 > + ^><°<>+^> 



1/ai - a 3 . 



(48) 



This is the Q-function for the nondegenerate parametric oscillator coupled to ordinary vacuum. On the 
other hand, in the absence of damping (7 = 0), Eq. (46) reduces to the form 



Q(a, a ,t) = exp 



— (tanh Kj t)(af3 + a* (3*) 



(49) 



which is the Q-function for the nondegenerate parametric amplifier. 

Next we obtain the Q-function for the single-mode generated by a degenerate parametric oscillator 
coupled to a single- mode squeezed vacuum reservoir from the Q-function for the NDPO (46). The Q- 
function for the single-mode can be expressed as 



a'b'c' + A'c' 2 + B'b 2 



Q(a,a*,t)= j d 2 /3Q(a,a*,(3,(3*,t), 
so that using Eq. (46) and the relation 

[d 2 aexp [- a'H 2 +6 / afc / a*+ 4^ 2 + R'n-* 2 ] = 1 n yp 

J(a' 2 -4A'B') 

the Q-function for the DPO coupled to a single-mode squeezed vacuum reservoir takes the form 

D 



-AA'B' 



Q(a, a*,t) 



exp 



a M 2 + ^-(a 2 + a* 2 



where 



bl-bj, 
1 



a = — 

V 

V 



(h + b 4 ) (h (h - 64) + 2b 2 b 3 ) - ft! (63 + b 3 ) 2 
(61 + 6 4 )(6 4 (bi - 64) + 2fe 2 6 3 ) + b±(b2 + h)' 



a' > (50) 



(51) 

(52a) 
(52b) 

(52c) 
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Upon integrating Eqs. (48) and (49) with respect to (3 by employing relation (50), one can also find 
the Q-function for the DPO in the absence of squeezed vacuum reservoir (r = 0) and in the absence of 
damping (7 = 0) to be 



and 



Q(a,a*,t) 



Q(a,a*,t) = 



n(ai + 03) 
sech 2 K7oi 



exp 



ai + a3 



exp [ - (sech 2 Kj t)(\a\ 2 )] , 



(53) 



respectively. 



5 Quadrature Squeezing 

In this section the intracavity quadrature fluctuations for the single-mode generated by the DPO as well 
as the signal-idler modes produced by the NDPO coupled to the two squeezed vacuum reservoirs using 
the pertinent Q-functions derived in the previous section are analysed. 

Here the first focus is the squeezing properties of the single- mode light. These properties could be 
described by two Hcrmitian operators defined as a\ — a? + a and 02 = i[a' — a). These quadrature 
operators obey the commutation relation [01,02] = 2i. The variance of these quadrature operators can 
be put in the form 



Aoi 



a l,2> 



(54) 



We now proceed to calculate the expectation values involved in expression (54). Applying the relation 



d 2 aQ(a, a* , t)A a (a, a*), 



(55) 



in which A a (a, a*) is the c-number equivalent of the operator A(a,a^) for the antinormal ordering, one 
arrives at 



/oo 
d 2 a Q(a, a* , t) a. 
-00 



Upon using the Q-function (51) for the single mode, the above equation can be expressed as 



■ exp 



aa a + — {a + a I 



ba 



b=0 



and on the basis of (50) for which d = and A' = B', one can verify that 

" CX P (2R' 2 



D d 
y/y db 



i 2 -A 2 ] 



Va 2 - A 2 



= 0. 



b=0 



In view of this result expression (55) reduces to 

(Aa l!2 ) 2 = 1 + 2(0+0) ± (o t2 ) ± (o 2 ). 



(56) 



Making use of the fact that the c-number equivalent of a* a for the antinormal ordering is a* a — 1 and 
applying relation (55) in evaluating all the expectation values in Eq. (56), we arrive at 



Aoi, 



aT A 



1. 



(57) 



Finally the quadrature fluctuations of the single-mode at any time t, in view of Eqs. (52), (46) and (45), 
take the form 



[Aoi 



Ai, 2 (e A5 *-l)+A 5 , A 3 , 4 (e- A6t -1) +A 6 



- 1. 
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At steady-state (t — > oo), the variances given above reduce to 



(\- \ 2 ^1,2 . ^3,4 

Aoi, 2 ) = — + — 
A5 Ae 



and with the aid of Eq. (26) one can rewrite these expressions as 

(Aai, 2 ) = 



Since for squeezed vacuum reservoirs 



1 



TV = sinh r, 

M = sinh r cosh r, 



(58) 



(59a) 
(59b) 



where r is the squeezing parameter taken to be real and positive for convenience, expression (59) takes 
the form 



(Aai, 2 )' 



Using (60) one can show that (Aai) 2 < 1, for 

r > In 

2 



1 



2^7o 



(60) 



(61) 



and (Aa 2 ) 2 > 1 for all r. This shows that the degenerate parametric oscillator coupled to a squeezed 
vacuum reservoir is in a squeezed state for the value of r specified by Eq. (61). 

In the absence of squeezing, i.e., r = 0, substitution of Eq. (26) into Eq. (57) leads to 



I _ «7o g— (7— 2«7p)t 

(Aar) 2 = (Aa 2 ) 2 = 1 



1 



' 2k7q \ 2 
v 7 / 



-4K7 n 



2^70 

7 



(1 



-4K7 n 



At steady-state and when the parametric oscillator is operating below threshold (7 > 2^70), this equation 
reduces to 



(AaO 2 = (Aa 2 ) 2 - 



1 



1_ (2^0)2 



in which both variances become greater than unity. Hence the single-mode in this case is not in a squeezed 
state. 

In the absence of damping (7 = 0), Eq. (57) reduces to 

(AaO 2 = (Aa 2 ) 2 =2n+l, 

where n = sinh 2 njot is the mean photon number for the single-mode. From these variances one can 
infer that the single-mode in this case is in a chaotic state as expected. Furthermore, in the absence of 
parametric interaction (re = 0), Eq. (57) could be expressed as 



(Aai, 2 ) 2 



2ai 2 -1 = 1- 



1 - e 



-7t 



T 2r 



1, 



and at steady-state these relations reduce to 

(Aa 1)2 ) 2 = e^ 2r , 



(62) 



which are the quadrature fluctuations of the squeezed vacuum reservoir A. 

Now we proceed to investigate the squeezing properties of the signal-idler modes produced by the 
NDPO coupled to the two squeezed vacuum reservoirs applying the Q-function (46). The squeezing 
properties of two-mode light can be described by two quadrature operators defined as 



Cl,2 



(63) 
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where 



&i = (o' + a), 
a 2 = i(cr — a), 



h = (6 f + b), 
b 2 = i(tf-b), 



(64a) 
(64b) 



and a (6) denotes the annihilation operator for the intracavity mode a (6). The quadrature operators c\ 
and c 2 satisfy the commutation relation [c\, c 2 ] = 2i. On account of these expressions, the variances can 
be expressed as 



(Ac li2 ) 2 = (c? i2 )-(c li2 r 



in which 



= ^1,2) + ^(R 2 ) + (ai,2, 01,2), 



(a,i,bi) = (aSi) - {ai)(bi}, 



(65) 



and i = 1,2. In particular, when a and 6 represent the signal and idler modes, respectively, it can be 
shown that 



(Aci, 2 ) 2 = i(A ai , 2 ) 2 + i(A 0l , 2 ) 2 + (ai, 2 , Si, 2 > 
= (Aai, 2 ) 2 + (01,2,61,2) 



(66) 



as (Adi j2 ) 2 = (A6i j2 ) 2 and (01,2) = (61.2) = (ci.a) = 0. In order to obtain the explicit form of Eq. (66), 
we proceed as follows. In view of expression (64) and (55), one can express that 



d 2 ad 2 f3 (a* + a)(/T + (3) Q(a,a* , /3* , (3,t). 



Then employing the Q-function (46) the above equation can be further expressed as 



(aih) = 



D f (a* + a) exp[ — bia* a + (a 2 + a* 2 ] 

J —00 7r ^ 



00 d 2 (3 



x / —(l3*+l3)ex P \-b 1 f3*f3+{b2a + b 3 a*)f3 
J-00 * 

+ (b 2 a*+b 3 a)0* + h 4 (p 2 +0* 2 )]. 



On setting K = 6 2 a + 63a*. 
(SiSi) = D 



— (a + a) exp — Oxce a + — {or + a ) 

7T L 2 



/ d d 
\dK + dK* 



00 12 



d 2 p 



-exp 



1, 



6i/3*/? + (6 2 a + & 3 a*)/3 + (& 2 a* + & 3 a)/3* + ^& 4 (/3 2 + /?* 2 ) 



so that performing the integration with respect to (3 on the basis of relation (50) and carrying out the 
differentiation we obtain 



(0161) = —(61 + 64) (62 + 63) 

from which it follows that 
D 



00 j2 r 

fl a / 2 1 *2 1 o * \ * 1 A r 2 1 *2\ 

(a + a + la a)exp — aa a -\ (a +a ) 

-no L 2 . 



(ai6i 



(6i + 6 4 )(6 2 + 6 3 ) 2 



00 ,2 



era 



^4 

aa*a + — (a 2 + a* 2 ) 



Next, integrating over a and carrying out the differentiation, we get 

D /2A + 2a\ 1 



(Mi) 



+ b 4 )(b 2 + 63) V a2 - A2 J Va 2 - A 2 ' 
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Making use of expression (47) along with Eq. (52) the above equation reduces to 

(ai&i) = ai - a 3 . 



A similar approach leads to 



Now Eq. (66) can be put as 



(02,62) = -(a 2 - 04). 



(Ac M ) 2 = 2a 1A - 1, 
and at this stage the variances are given by 



K70 ± 7(iV =f M + 1) 



2k7q ± 7 



(7 - e T(2K7o±7)*j + e T(2K7o±7)< 



1. 



Finally the quadrature fluctuations of the signal- idler modes at any time t, in view of Eq. (59), take the 
form 



(Ac x ) =1- 



[Ac 2 ) =1- 



1 _ g-(7+2K7o)* 

2 _ g-(7-2K7o)* 



7 e 



1 - 



7 + 2k7 

7e +2r 



7 - 2^70 



< 1, 



> 1. 



(67a) 
(67b) 



Hence the signal-idler modes generated by the NDPO coupled to the USVR, when operating below 
threshold (7 — 2«7 > 0), are in squeezed states for all values of r. 
At steady-state (t — > 00), Eq. (67) can be put in the form 



Acr 



(Ac 2 ) = 



7 



7 + 2«7o 



e~ 2r < 1, 



7 



7 - 2k7q 



,+2r 



> 1. 



(68a) 
(68b) 



This equation clearly shows the possibility of a very large amount of squeezing (approaching 100%) 
below the standard quantum limit in the one quadrature at the expense of enhanced fluctuations in the 
other quadrature, where in this case the standard quantum limit is taken to be */ (Aci) 2 y (Ac 2 ) 2 = 1. 
In addition, at threshold (7 = 2^70), one obtains 



(Acr) 2 
(Ac 2 ) 2 



(69a) 
(69b) 



In the absence of squeezed vacuum reservoirs (r = 0), expression (68) becomes 



(Ac li2 ) 2 



7 ± (2^0 e 



(7±2«7o)t 



7 ± 2k7q 



SSI- 



(70) 



This shows that the signal-idler modes produced by the nondegencratc parametric oscillator in the absence 
of squeezed vacuum reservoirs are also in squeezed states. At steady-state and at threshold, these relations 
reduce to 



(Acr) 2 
(Ac 2 ) 2 



1 

2' 
00. 



(71) 



In this case one can easily see that there is only a 50% reduction of noise below the vacuum level. By 
comparing Eqs. (69) and (71) we can conclude that coupling of the NDPO to the squeezed vacuum 
reservoirs is essential for the generation of a larger amount of squeezing. 
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In the absence of damping (7 = 0), expression (67) reduces to 



Ac 2 j2 = e T2K70t $ 1, (72) 

which are the quadrature fluctuations of the signal-idler modes produced by the nondegenerate paramet- 
ric amplifier. This indicates that the nondegenerate parametric amplifier coupled to ordinary vacuum 
reservoirs also generates squeezed states. 

Finally, when there is no parametric interaction inside the cavity (k = 0), Eq. (67) takes the form 

(Ac li2 ) 2 = 1 - (1 - e"*) [l T e* 2r ] $ 1, (73) 

which, at steady-state, leads to 

(Acr, 2 ) 2 = e ±2r , (74) 

which are the quadrature fluctuations of the reservoir modes A and B. Upon comparing the relations (71) 
and (74) with (69), one can see that the quadrature variances at steady-state and at threshold are the 
product of the variances of the NDPO coupled to ordinary vacuum and the variances pertaining to the 
squeezed vacuum reservoirs. Furthermore upon comparing expressions (62) and (74) one can observe 
that at steady-state the variances of a signal mode squeezed vacuum reservoir as well as those of two 
independent squeezed vacuum reservoirs are the same. 



6 Conclusion 

We have derived the master equation for the signal-idler modes produced by the nondegenerate parametric 
oscillator coupled to two uncorrclatcd squeezed vacuum reservoirs and consequently the Fokker-Planck 
equation. We have solved the pertinent Fokker-Planck equation which is a second order differential 
equation applying the propagator method [28] and obtained a compact form of the Q-function of the 
optical system coupled to two independent squeezed vacuum reservoirs. We have also deduced the Q- 
functions for a NDPO coupled to ordinary vacuum reservoirs, degenerate parametric oscillators coupled to 
a squeezed vacuum reservoir and an ordinary vacuum reservoir, and for the nondegenerate and degenerate 
parametric amplifiers from the Q-function for the NDPO coupled to the two USVR. 

In general the Q-function can be used to evaluate the expectation values of antinormally ordered 
operators as well as photon number distributions for the NDPO and other similar optical systems. 

We have calculated the nonlinear quantum quadrature fluctuations of the signal-idler modes gen- 
erated by a nondegenerate parametric oscillator below threshold coupled to two uncorrclatcd squeezed 
vacuum reservoirs, using the Q-function. Although it is a well known fact that quantum noise can not be 
eliminated, we have shown that the signal-idler modes produced by the optical system are in a two-mode 
squeezed state at any time t. More interestingly, we have shown that at steady-state and below threshold 
it is possible to generate an optimal squeezing in one of the quadratures below the standard quantum 
limit at the expense of enhanced fluctuations in the other quadrature so that the Hciscnberg uncertainty 
principle remains valid. Furthermore calculation of the quadrature fluctuations at threshold clearly shows 
that it is possible to produce an arbitrarily large squeezing (approaching 100%) in one of the quadratures 
with an infinitely large noise in the other quadrature. We have also shown that the degenerate parametric 
oscillator could be in a squeezed state for a squeezing parameter above a certain value when it is coupled 
to a squeezed vacuum reservoir. 

We have shown that the coupling of the optical system to the squeezed vacuum reservoirs is essential 
in order to get a more suppressed noise in one of the quadratures. 

Finally we have calculated the quadrature fluctuations for the nondegenerate parametric amplifier 
coupled to ordinary vacuum reservoirs and verified that it also generates squeezed states. 
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